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The medium temperature T dependence of jet transport coefficient ĝ is studied via nuclear modification factor 
Raa(pr) and elliptical flow parameter vo(pr) for large transverse momentum pr hadrons in high-energy 
nucleus-nucleus collisions. Within a next-to-leading order perturbative QCD parton model for hard scatterings 
with modified fragmentation functions due to jet quenching controlled by ĝ, we check the suppression as well 
as the azimuthal anisotropy for large pr hadrons, and extract ĝ by global fits to Raa(pr) and v2(pr) data 
in A + A collisions at RHIC and the LHC, respectively. Numerical results from the best fits show that ¢/ 7 
goes down with the local medium temperature T in the parton jet trajectory. Compared with the case of a 
constant ĝ/ T°, the going-down T dependence of ĝ / T? makes a hard parton jet to lose more energy near T. and 
therefore strengthens the azimuthal anisotropy for large pr hadrons. As a result, v2(pr) for large pr hadrons 
is enhanced by about 10% to fit data better at RHIC/LHC. Considering the first-order phase transition from 
QGP to the hadron phase and additional energy loss in the hadron phase, v2(pr) is again enhanced by 5-10% at 


RHIC/LHC. 
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I. INTRODUCTION 


The suppression and azimuthal anisotropy of the high 
transverse momentum (pr) hadrons are two valuable evi- 
dences for the existence of the quark-gluon plasma (QGP) 
that might be created in high-energy nucleus-nucleus colli- 
sions performed at both the Relativistic Heavy-Ion Collider 
(RHIC) [1-3] and the Large Hadron Collider (LHC) [4— 
12]. When high energy partons propagate through the color- 
deconfined QGP medium, they will encounter multiple scat- 
terings and lose their energy via the medium-induced gluon 
radiations. Then the final state hadrons observed in nucleus- 
nucleus (A + A) collisions will be suppressed compared to 
those observed in proton-proton (p + p) collisions. Gener- 
ally, the suppression strength is given by the nuclear modi- 
fication factor R4,4(pr) defined as a ratio of single hadron 
spectrum in A + A collisions over that in p + p collisions. 
In a typical non-central A + A collision, the initial geometric 
anisotropy can be converted into the azimuthal anisotropy in 
the gluon density distribution of the produced QGP medium, 
which leads to the azimuthal anisotropy of the total energy 
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loss for energetic jets due to the path length and the gluon 
density dependence of the jet energy loss. To characterize 
this anisotropy, one can introduce the elliptic flow coeffi- 
cient v2(pr) which is defined as the second order Fourier co- 
efficient in the azimuthal angular distribution of final state 
high pr hadrons. Both of the two observables R4a(pr) 
and v2(pr) for large pr hadrons are the consequence of jet 
quenching or energy loss [13—20], which are expected to give 
a consistent jet quenching description. 

The strength of jet energy loss is controlled by jet transport 
coefficient ĝ, which is proportional to the medium gluon num- 
ber density p and defined as the average transverse momen- 
tum broadening qr squared per unit length for jet propagating 
inside the medium [21], 


î= p | ahh (1) 


Quantitative extraction for the energy loss parameter was first 
performed by JET Collaboration utilizing different theoretical 
models and different approximations compared to experimen- 
tal data for single hadron productions at RHIC and the LHC 
[22]. For simplicity, ĝ/T? is generally assumed as a con- 
stant to study bulk matter evolution [23] and the suppression 
of large pr single hadron and dihadron production [24, 25]. 
The jet transport coefficient and the mean free path at the ini- 
tial time are simultaneously extracted for jet energy loss [26]. 
Comparing the extracted ¢/T? for the different initial tem- 
peratures in the center of QGP between RHIC and the LHC 
cases, numerical results hint that ¢/T? decreases slightly with 
the increasing of the medium temperature [22—26]. MAR- 
TINI [27], MCGILL-AMY [28] and other theoretical stud- 
ies [29, 30] also unfold similar conclusions. In fact, pertur- 


bative studies with resummed hard-thermal loops in finite- 
temperature QCD give rise to an additional temperature de- 
pendence of ¢/T°? for a fixed value of the strong coupling con- 
stant [31, 32]. CUJET model [33, 34] considers there might 
exist a strong dependence of å/T? on temperature and try to 
give a systemic description on R4,4(pr) and vo(pr) simul- 
taneously with a Gaussian-like temperature dependence form 
of G/T? [35, 36] within the opacity expansion energy loss for- 
malism [37]. There are also many other descriptions and de- 
velopments for the jet quenching parameter, such as the radia- 
tive corrections to ĝ [38-40] and nonperturbative calculations 
for it using the AdS/CFT correspondence at strong coupling 
in string theory [41-43] and lattice approaches [44—47]. Re- 
cently, with newly-developed Bayesian analysis, JETSCAPE 
has studied the medium temperature, virtuality, and jet energy 
dependence of ĝ via single hadron suppression at RHIC and 
the LHC energies [48]. Meanwhile, the LIDO model [49] and 
JETSCAPE [50] also extracted the ĝ value with two types of 
observables, single inclusive hadron and jet suppression. In 
Ref. [51], with non-parametric prior distribution of ĝ, using 
single hadron production, dihadron and y-hadron correlation 
data calibrated the temperature-dependent g. All studies point 
out that the g/T? should have a larger value at critical tem- 
perature Te. 


In this paper, we study on the additional temperature de- 
pendence of ¢/T? via comparing theoretical calculations with 
experimental data for both Raa(pr) and v2(pr) at large pr 
at RHIC and the LHC. To reveal a clear tendency for the 
additional temperature dependence, we assume a linear or 
Gaussian-distribution form for the temperature dependence 
of g/T? within a high-twist energy loss formalism [52-54]. 
A (3+1)d ideal hydrodynamic description of the bulk matter 
evolution is used for the medium expansion, which is out- 
putted in references [55, 56] for Au + Au collisions at 200 
GeV and Pb+ Pb collisions at 2.76 TeV. The initial condition 
for the ideal hydrodynamic equations is fixed so that the fi- 
nal bulk hadron spectra from the experiments are reproduced. 
In order to fit both R4a(pr) and v2(pr) simultaneously, we 
first consider only the QGP phase for jet energy loss and then 
the hadron phase contribution [23] is also included. Our cal- 
culations give a good description to R4a,(pr) for different 
temperature dependence schemes of ¢/T?, while theoretical 
results for v2(pr) underestimate the experimental data. How- 
ever, comparing to the case with a constant ¢/T°, we find that 
the going-down temperature dependence of @/T® gives about 
10% rise to v2(pr) in the QGP phase, and an additional 10% 
rise at RHIC and 5% rise at the LHC when hadron phase con- 
tribution is included. 


The rest of the paper is organized as follows. We first 
review the next-to-leading-order (NLO) perturbative QCD 
(pQCD) parton model with the medium modified fragmen- 
tation functions in Sec. II. Then shown in Sec. III and Sec. IV 
are our numerical results fitting to R44(pr) and v2(pr) data 
for linear and Gaussian temperature dependence of ¢/T® in 
QGP phase, respectively. In Sec. V the hadron phase contri- 
bution is also included for the linear temperature dependence 
of G/T°?. Finally, we conclude in Sec. VI with a summary. 


Il. NLO PQCD PARTON MODEL WITH MODIFIED 
FRAGMENTATION FUNCTIONS 


Within the NLO pQCD parton model, the collinear fac- 
torized differential cross-section of single hadron production 
in p + p collisions can be factorized into the convolution of 
parton distribution functions (PDFs), short-distance partonic 
cross-sections and fragmentation functions (FFs) [57, 58], 
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where fajp(£a, H?) is the parton distribution function for the 
parton a with the momentum fraction x, from a free nucleon, 
and we will use CT14 parametrization [59]. The fragmenta- 
tion function D} (ze, u?) for a parton in vacuum is given by 
the AKK parametrization [60], in which z, is the momentum 
fraction carried by the outgoing hadrons from the parent par- 
ton c. da(ab — cd)/dé is the parton-parton hard scattering 
cross section at LO a2. In our numerical simulations, the par- 
tonic scattering cross sections will be computed up to NLO 
implied in O(a). The NLO corrections include 1-loop con- 
tributions to 2 — 2 tree level and 2 — 3 tree level contribu- 
tions. More detailed discussions on NLO calculations could 
be found in the reference [61]. 

In A + A collisions, the cross-section for single hadron 
production at high transverse momentum is given by [62, 63], 
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where t4(7) = f pa(7)dz is the nuclear thickness function 
given by the Woods-Saxon distribution and is normalized as 
[@rta(r) = A. fasa(ta, H’, T) is the nucleus-modified 
parton distribution function that is assumed to be factorized 
into the parton distributions in a free nucleon fa/y (£a, yl?) 
and the nuclear shadowing factor S,./4(a, 7, r) [64, 65], 


Z 
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where Z is the proton number of the nucleus and A is the 
nucleus mass number. Assuming that the shadowing is pro- 
portional to the local nuclear density, the shadowing fac- 
tor SajAlTa, u?,7) can be obtained using following form 
[66, 67], 


1] Ata (r) 
J Prita)? 
where SajA(Za, u?) is given from the EPPS16 [68]. 


SajA (Ea, WT) = 1 +F [Seat ) (5) 


The medium-modified fragmentation functions ri, can be 
calcualted as follows [62, 63, 69-71]: 


Digar AEs) = (1— e7'%2)) E Ze DACA, p?) 
Zh - 
i (0) DB 1?) +N) Di (ze u’), © 


where z. = pr/pre is the momentum fraction for a parton 
fragmenting into a hadron in vacuum. zc’ = pr/(pre—-AE-) 
is the rescaled momentum fraction, and denotes that a parton 
with pr, propagating through the medium loses energy AE, 
and fragments into a hadron with pr. zg’ = pr/(AE./(N,)) 
is the momentum fraction for a radiated gluon fragmenting 
into a hadron. (N,) is the radiated gluon number. 

The parton energy loss caused by medium-induced gluon 
radiation can be calculated in the higher-twist (HT) approach 
[52-54]. For a light quark c with initial energy F, the radia- 
tive energy loss AE, can be calculated as, 
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where C4 = 3, a, is the strong coupling constant, and Ip is 
the transverse momentum of the radiated gluon. We assume 
the energy loss of a gluon is simply 9/4 times that of a quark 
due to the different color factors for the quark-gluon vertex 
and gluon-gluon vertex [52]. The average number of radiated 
gluons from the propagating hard parton is calculated as [72], 
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The HT formalism contains the transverse momentum lr of 
the radiated gluon, which also indicates changes in the trans- 
verse momenta of the partons [15, 73, 74]. In our numerical 
simulations, we adopt the small angle approximation within 
the collinear factorization theorem according to Eq. (3). Con- 
sequently, we focus solely on the effect of energy loss and 
assume the parton direction remains unchanged in the frag- 
mentation functions. Such an approximation has been used in 
many current jet energy loss formalisms and has successfully 
explained experimental data [49, 50, 75-78]. 

The parton energy loss and the number of the radiated glu- 
ons are both controlled by the jet transport parameter ĝ [21]. 
According to Eq. (1) for ¢ proportional to the medium gluon 
density p, one can simply assume a constant value of the 
scaled jet transport parameter [25, 26], 
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where p” is the four momentum of the parton, u” is the local 
four flow velocity of the fluid, T' is the local temperature of 
the medium and Tọ is a reference temperature taken as the 


highest temperature at the center of the medium at the initial 
time To. 

For an additional T dependence of å/T?, one can simply 
assume a linear form like G/T? ~ aT + b. In the following 
actual calculations, we write the linear form as, 
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T= T? 
We also check the additional T dependence of ĝâ/T° by Gaus- 
sian form, 
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The parameters ĝo, ĝe, and or are introduced to adjust the 
strength of the additional temperature dependence. Te = 170 
MeV is the critical temperature. When T = To for Eq. (10), 
and 07, = oo for Eq. (11), the two equations both go back to 
Eq. (9), respectively. 

To describe jet quenching in high-energy nucleus-nucleus 
collisions, one needs to provide the space-time evolution of 
the jet transport coefficient in Eq. (9, 10, 11) along the par- 
ton propagation. In our studies, the dynamical evolution of 
the medium that governs the space-time evolution of the lo- 
cal temperature T and flow velocity u is obtained using the 
(3+1)-dimensional hydrodynamic model [55, 56]. This model 
gives results on the transverse dynamics of the bulk medium 
in A+ A collisions with the initial conditions. Furthermore, 
the model includes the first order phase transition between the 
QGP phase and the hadron phase at Te = 170 MeV and gives 
the hadron phase fraction f(r) which is defined as, 


0 if T>170 MeV, 
f(® =< O~1 if T=170 MeV, (12) 
1 if T<170 MeV, 


where 7 is the jet local position. As the geometry position 
moves closer to the periphery of the medium or the medium 
evolution time becomes longer, the hadronic phase fraction 
gradually increases from 0 to 1 at Te = 170 MeV, determined 
by the proportion of hadron and parton number density [55, 
79, 80]. The energy loss of the jet propagating through both 
the QGP phase and the hadronic phase can be described by 
the higher-twist approach simultaneously, except that the ĝ in 
the separate phase is different. To include the contributions to 
q from both QGP and hadron phases, we change Eq. (9), (10) 
and (11) to be 


dh 
Dral (13) 
where ĝp is the jet transport parameter for the hadronic phase. 
By combining Eq. (12), one can get that for studies exclu- 
sively concerning the QGP phase, we consider the pure par- 
tonic medium at T > 170 MeV, along with the QGP fraction 
(1 — f) in the mixed phase at Te = 170 MeV. For studies 
on the hadronic phase, we only need to consider the By f 
term that accounts for the contributions from the hadronic 
phase during the mixed phase, as well as the entire hadronic 
medium when the temperature is below 170 MeV, until the 


system reaches dynamic freeze-out. When the hadron phase 
is also considered for jet energy loss, the extracted jet trans- 
port parameter for the QGP phase is reduced due to a long 
evolution time of the mixed phase, as shown in Ref. [23]. 

The jet transport parameter in the hadron phase can be writ- 
ten as [23], 


_ an 2 
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where Gy ~ 0.02 GeV7/fm is the extracted jet transport pa- 
rameter at the center for the cold nucleonic matter of a large 
nucleus, and py ~ 0.17 fm~® is the nucleon density at the 
center of the large nucleus [52]. pm and ppg are the meson 
and baryon density in the hadronic resonance gas at a given 
temperature, respectively. Factor 2/3 represents the ratio of 
the constituent quark number of the meson and baryon. The 
hadron density at a given temperature T’ and zero chemical 
potential is expressed as [23], 
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where ņa = + for meson (M)/baryon (B). In the fol- 

lowing calculations, hadron resonances with mass 

below 1 GeV are included: n kinds of meson, 


nt, n,n, K+, K~, K®, K}, n,n, pt, p~, p, K**, 
K*-, K*, K*, w; and 2 kinds of baryon, p,n. Here, we 
ignore the cöntribütion of the anti-nucleons to ĝa, which is 
less than 3%. 


II. LINEAR TEMPERATURE DEPENDENCE OF ĝ/T? IN 
QGP PHASE 


With the spectrum in p + p collisions as a baseline, the 
nuclear suppression factor R4a(pr) for single hadron pro- 
duction in A + A collisions can be expressed as [71, 81], 


dN} , /dyd*py 
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where T44 (b) = f d?rt.a(F)tg(?+6) is the overlap function 
of two tite nuclei. 

The anisotropy of the final state hadrons in transverse mo- 
mentum can be quantified by the Fourier expansion of the 
hadrons distribution in the azimuthal angle. We focus on the 
second Fourier coefficient, namely elliptic anisotropy coeffi- 
cient v2(pr), which can be written as [82-87], 


i, docos(2¢)dN} , /dyd?prdo 
i dod N}, , /dyd?p rdo Í 


where ¢ is the jet azimuthal angle between the jet’s propagat- 
ing direction and the impact parameter. 

In this section, we will use Eq. (10) and (13) with ĝa = 0 to 
consider the linear temperature dependence of ¢/T? in QGP 
phase. x? fitting to both R44(pr) and v2(pr) for hadrons 


v2(pr) = (17) 


in the middle rapidity region will be performed for different 
introduced parameters, which is given by, 


=i (Vin = Vexp ay ee + Orta) y (18) 
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where Vin and Vexp denote the theoretical and experimen- 
tal results, and osys and Ostat give systematic and statisti- 
cal errors for the data, respectively. For a global fit for both 
Raa(pr) and v2(pr), the data number N for the degree of 
freedom (d.o.f.) is the sum of Raa(pr) and v2(pr) data 
numbers. The y?/d.o.f value is minimized to near unity to 
find the best-fit temperature dependence of @/T® [36]. In the 
x?/d.o.f calculations, we only select the experimental data 
points with pr > 7.5 GeV/c at both RHIC and the LHC en- 
ergies to ensure the validity of pQCD parton model. 


A. Fit R44 and v2 at RHIC 


Current researches all pointed that the @./T? at critical 
temperature has a larger value, and the value of go /Tj at 
the highest temperature is smaller. Therefore we first choose 
Gc/T3 € [3.0,9.0] and ĝo/TË € [0.2,5.6] with bin size 0.3 
and get 399 couples of (G./T?, Go/T) for Eq. (10) and (13) 
with ĝa = 0. To find the limit value of go/T3, we could 
expand it to 0 phenomenally. So we make 420 times of cal- 
culations for the suppression factor R44(pr) as a function 
of pr for single hadrons produced in the most central 0 - 5% 
Au + Au collisions at \/snn = 200 GeV. Each result for 
Raa(pr) with a given couple of (¢./T?, Go/Tj) provides 
one value of x? /d.o. f to fit to experimental data of Raa(pr) 
[1, 2]. As shown in Fig. 1 (a) is such a 2-dimensional figure 
for x*/d.o.f as a function of (¢./T3, Go/T}). The different 
colors represent different fitting values. The x? /d.o.f value 
is minimized to near unity to find the best-fitting couples of 
(G-/T?, ĝo/T). 

The best fittings given by the red region show that the sin- 
gle hadron R44(pr) is more sensitive to the value of ¢./T? 
than ĝo/T9. But the fitting fails to get one unique couple of 
(åe/T3 , ĝo /T9) for an explicit dependence form of g/T* on 
T only through the constraint of single hadron R4 (pr). To 
show the different linear temperature dependence of ¢/T® for 
the same suppression of single hadrons, we draw in Fig. 1 (b) 
the gray dashed curves for @/T® as a function of T, which are 
constrained by the best-fitting region of x? /d.o.f. 

Among the gray dashed curves, we choose one horizon- 
tal line (blue) for a constant ¢/T? with (g./T?, Go/T3) = 
(5.6, 5.6) and one ae line (red) for a linear T' dependence 
of ¢/T? with (G-/T?, Go/T§) = (7.2, 0.0). Using these two 
couples of (@./T3, ĝo /TÈ ), we get almost the same R4a(pr) 
in Fig. 1 (c) which shows that the single hadron suppression 
is a consequence of total jet energy loss and is not sensitive to 
the T dependence of g/T® in central Au + Au collisions. 

Due to the dependence on the jet path length and the 
medium density in the jet trajectory inside the hot medium, jet 
energy loss in noncentral Au + Au collisions has azimuthal 
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Fig. 1. (Linear-dependence) Panel (a): The x?/d.o. f analyses for single hadron R44 (pr) as a function of ¢./T? and go/T@ from fitting to 
PHENIX data [1, 2] in the most central 0 - 5% Au + Au collisions at ,/snn = 200 GeV. Panel (b): The scaled dimensionless jet transport 
parameters ĝ/ T? as a function of medium temperature T from the best fitting region of the panel (a). Panel (c): The single hadron suppression 


factors R4a(pr) with couples of (ĝe /T3, ĝo /TÈ) = 
[1, 2] data. 


(7.2, 0.0) (red solid curve) and (5.6, 5.6) (blue dashed curve) compared with PHENIX 
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Fig. 2. (Linear-dependence) The x° /d.o.f analyses for single hadron R.a(pr) (panel (a)) and elliptic flow v2(pr) (panel (b)) as a function 
of G/T? and ĝo/TÈ from fitting to experimental data [1-3] in 20 - 30% Au + Au collisions at \/snn = 200 GeV. The global x? /d.o. f 


fitting results for both R.44(pr) and v2(pr) are shown in the panel (c). 


anisotropy. The hadron suppression depends on the azimuthal 
angle with respect to the reaction plane, thus leading to az- 
imuthal anisotropy in high pr hadron spectra. The same en- 
ergy loss mechanism permits to make a global fit to constrain 
(âe / T2, ĝo /TÈ) with both the suppression factor R4(pr) 
and the elliptic flow parameter vo(pr) for large pr hadrons 
in noncentral A + A collisions. With the same couples of 
(G-/T?, Go/Te) as in 0 - 5% centrality, we simultaneously 
make 420 times of calculations for Ra ,a(pr) and v2(pr) as a 
function of pr to fit to the experimental data [1-3] in 20 - 30% 
Au-+ Au collisions, and get the x? /d.o. f results for R44 (pr) 
in Fig. 2 (a) and vo(pr) in Fig. 2 (b), respectively. The 
x?/d.o.f fitting for Raa in noncentral collisions is similar 
to that in central collisions. This means that only RAA (pr) 
constraint does not give an explicit dependence form of ¢/T? 
on T. The y?/d.o.f fitting for v2 in Fig. 2 (b) shows that the 
data of elliptic flow v2(pr) favor larger ĝe /T3 and are almost 
insensitive to Go/T@. A global x?/d.o.f fit is performed for 
both R44(pr) and v2(pr) in Fig. 2 (c) in which the limited 


yellow region is found to constrain (G./T?, Go/T@). 


Shown in Fig. 3 (a) is the scaled dimensionless jet trans- 
port parameters ¢/T®? as a function of medium temperature 
T from the best fitting region of global y? fits of Fig. 2 
(c). The blue dashed curve is also for the constant de- 
pendence case with (g./T?,Go/Tj) = (5.6,5.6), and the 
red solid curve for a linear T dependence of ¢/T? with 
(G./T?, Go/Tj) = (6.9,0.0). These two dependence forms 
almost give the same R44(pr) as shown in Fig. 3 (b), which 
is similar to the situation in central collisions. However, these 
two different dependence of ¢/T® gives different contribu- 
tions to v2(pr) as shown in Fig. 3 (c). Numerical results 
show that the linearly-decreasing T dependence of ¢/T® with 
(G-/T?, âo /T8) = (6.9, 0.0) makes an enhancement by 10% 
for ve(pr) comparing to the constant dependence case with 
(G./T?,Go/T3) = (5.6,5.6). Such a linearly-decreasing T 
dependence of ĝâ/T? with (¢./T?, Go/Tj) = (6.9, 0.0) means 
more energy loss taking place near the critical temperature Te. 
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Fig. 3. (Linear-dependence) Panel (a): the scaled dimensionless jet transport parameters ĝâ/T? as a function of medium temperature T from 
the best fitting region of global x° fits of Fig. 2 (c) in 20 - 30% Au + Au collisions at \/snn = 200 GeV. The single hadron suppression 


factors R.44(pr) and elliptic flow v2(pr) are shown in panel (b) and (c), respectively, with couples of (q-/T?, Go/Té) = 
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solid curve) and (5.6, 5.6) (blue dashed curve) compared with PHENIX [1-3] data. 
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data. 


B. Fit Ry, and v2 at the LHC 


Similarly, we present our relevant results for Pb + Pb col- 
lisions at \/snN = 2.76 TeV. Here we choose ĝ./T? € 
[1.2,8.1] and go/T? € [0.2,6.5] with the same bin size 0.3 
and further include the ĝo/T = 0.0, thus get 552 couples 
of (Gc/T?, Go/T 3) for Eq. (10) and (13) with ân = 0. The 
x?/d.o.f results for central 0 - 5% Pb + Pb collisions are 
performed on single hadron suppression factors, as shown in 
Fig. 4 (a). The best-fitting contour is similar to that in Fig. 
1 (a) but with a smaller ¢./ T: With a constant form and a 
linear form for ĝâ/T?, we again obtain the same single hadron 
suppression as shown in Fig. 4 (b). 


For 20 - 30% Pb + Pb collisions, the y?/d.o.f fitting for 
only Raa(pr) or ve(pr), and the global fitting for both of 
them are shown in Fig. 5 (a), (b) and (c), respectively. Sim- 
ilar to noncentral Au + Au collisions, both the separated 
x?/d.o.f fitting for R44(pr) and vo(pr) can not give clear 
constrain on the T dependence of ĝ/ T. However, the differ- 


ence from \?/d.0.f(Raa), x7/d.0.f(v2) shows that data of 
v2(pr) prefer larger jet energy loss near T, and are insensitive 
to the change of ĝo /Tj. Consequently, the global fits for both 
Raa(pr) and v2(pr) make the constrain to some extent on 
the T dependence of @/T® as shown in Fig. 5 (c), similarly to 
Fig. 2 (c). 


Choosing y?/d.o.f < 1.6 in Fig. 5 (c) for the best fitting, 
one can get the curves for the T dependence of G/T? in Fig. 
6 (a). We again get a tendency for g/T® to decrease with the 
increasing local temperature in the jet trajectory. Among the 
best-fitting values, selecting (G-/T?, Go/T$) = (5.4, 0.0) for 
the linear T dependence and (4.1,4.1) for the constant de- 
pendence, we calculate the R44(pr) and v2(pr) as a func- 
tion of pr shown in Fig. 6 (b) and (c), respectively. Two 
almost same R.44(pr) are obtained, while the va(pr) has an 
enhancement by 10% due to the larger jet energy loss near 
T, for the linearly-decreasing T dependence of ¢/T® at the 
LHC. 


No matter whether in Pb + Pb collisions or Au + Au colli- 
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fitting results for both R44 (pr) and v2 (pr) are shown in the panel (c). 
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Fig. 6. (Linear-dependence) In 20 - 30% Pb + Pb collisions at \/snn = 2.76 Tey, the scaled dimensionless jet transport parameters 
â/T? as a function of medium temperature T from the best fitting region of global x° fits of Fig. 5 (c) are shown in the panel (a). The 
single hadron suppression factors Raa(pr) and elliptic flow v2(pr) are shown in the (b) and (c) panels, respectively, with couples of 
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sions, R4a(pr) and v2(pr) are both more sensitive to the jet 
energy loss near critical temperature T, than near the initial 
highest temperature To. Relatively, the data of v2(pr) prefer 
a larger value of g./T. Furthermore, the anisotropy of final 
state hadrons at high transverse momentum can be strength- 
ened up to 10% by increasing the jet energy loss near Te with 
the linearly-decreasing T dependence of ¢/T?. 


C. Jet energy loss distribution 


Given a parton jet with any creation site and any moving 
direction in the initial hard scattering, let us see the jet energy 
loss distribution when propagating through the hot medium. 
The average energy loss rate in the jet trajectory is given by, 


_ fideo f @rta(*)ta(F + b)dAE(F + 


( = nr) /dr 
dr J do f @rta(*tp(F+ 6) 


? 


(19) 


(5.4, 0.0) (red soild curve) and (4.1, 4.1) (blue dashed curve) compared with experimental data [4, 5, 9, 10]. 


where AF(7r) is given by Eq. (7), 7 is the initial creation 
point for an energy-given jet, and ři is the unit vector of the 
jet moving direction @ which is the same as that in Eq. (17). 
The average accumulative energy loss for the jet traversing 
through the medium is then given as, 


J do | Prta(F)ta(F +b) [7 dr SA 
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Shown in Fig. 7 are the average accumulative (solid 
curves) and differential (dashed curves) energy loss for one 
10 GeV jet with (G/T ?, Go/T3) = (6.9, 0.0) (red curves) and 
5.6, 5.6) (blue curves) in 20 - 30% Au + Au collisions at 
/8nn = 200 GeV (panel (a) and (b)), and for one 100 
GeV jet with (g./T?, Go/Tj) = (5.4,0.0) (red curves) and 
4.1,4.1) (blue curves) in 20 - 30% Pb + Pb collisions at 
V/S8NN = 2.76 TeV (panel (d) and (e)), respectively. The 
medium temperature as a function of time at the center point 
x,y) = (0,0) for the two collision systems are shown in the 
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Fig. 7. (Linear-dependence) Panels (a) and (d): the average accumulative energy loss for one 10 GeV jet in 20 - 30% Au + Au collisions 
at \/snn = 200 GeV and for one 100 GeV jet in 20 - 30% Pb + Pb collisions at \/snn = 2.76 TeV, respectively. Panels (b) and (e): the 
corresponding differential energy loss. Red solid or dashed curves are for the linearly-decreasing T dependence of (q/ T°, and blue for the 
constant (G/T°. Panels (c) and (f): the medium temperature as a function of time at the center point (x, y) = (0,0) for the two collision 


systems, respectively. 


lower panels (c) and (f), respectively. 

When the jet passes out of the critical region from QGP 
to the hadron phase, it is over for the jet to accumulate the 
lost energy, as shown in Eq. (13) with q@, = 0. For the well- 
chosen two cases of the constant dependence and the linearly- 
decreasing T dependence of ¢/T°, the final total energy loss 
is exactly similar, as shown in Fig. 7 (a) and (d). In the mean- 
time, the peak of the jet energy loss distribution (dA E/drT) 
along the jet path is “pushed” to move to critical tempera- 
ture T, nearby due to the linearly-decreasing T dependence 
of g/T® (red dashed curves) compared to the constant depen- 
dence (blue dashed curves), as shown in Fig. 7 (b) and (e). 
More energy loss happens on the approaching critical tem- 
perature and enhances the final hadron azimuthal anisotropy, 
so v2(pr) is strengthened up for the linearly-decreasing T 
dependence case, as shown in Fig. 3 (c) and Fig. 6 (c). 

In order to clearly illustrate the enhanced azimuthal 
anisotropy, we define the energy loss asymmetry below, 


(AE(1)) 9°"? — (AE(1)) 9 

(AE(r))e="/2 + (AE(r))e=0" 
where (AE(r, ¢)) is given by Eq. (20) in which the ¢ in- 
tegration for the azimuthal average is removed. On average, 
a parton jet will encounter the most energy loss due to the 
longest path length at ¢ = 7/2 and the least at 6 = 0. Shown 
in Fig. 8 (a) and (b) are such energy loss asymmetries for an 


A(AE(r)) = (21) 


energy-given parton jet in 20 - 30% Au + Au at 200 GeV and 
Pb + Pb collisions at 2.76 TeV, respectively. The red solid 
curves are for the case of the linearly-decreasing T depen- 
dence of g/ T®, while the blue solid curves are for the constant 
case. The former is 10% larger than the latter in both panels, 
which is similar to the enhancement for v2(pr) in Fig. 3 (c) 
and 6 (c). 

Due to the medium temperature evolution, different T de- 
pendence of jet transport coefficient gives different energy 
loss distribution for jet propagation. The large pr hadron sup- 
pression R44 is a consequence of total energy loss and inde- 
pendent of the jet energy loss distribution. However, com- 
pared to the constant case for a given total energy loss, the 
linearly-decreasing T dependence of ¢/T® makes an energy 
loss to redistribute and leads to more energy loss happening 
near the critical temperature and, therefore, more stronger az- 
imuthal anisotropy for hadron production. 


IV. GAUSSIAN TEMPERATURE DEPENDENCE OF ¢/T® 
IN QGP PHASE 


In the last section, numerical results for the linear T de- 
pendence assumption show that ĝâĝ/T goes down with the 
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Fig. 8. (Linear-dependence) The energy loss asymmetry between the jet propagating direction of ¢ = 7/2 and ¢ = 0 for one 10 GeV jet with 
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TeV (panel (b)). 
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Fig. 9. (Gaussian-dependence) Panel (a): The x” /d.o.f analyses for single hadron Raa(pr) as a function of o7:/T? and go /T§ from fitting 
to experimental data [1, 2] in the most central 0 - 5% Au + Au collisions at snn = 200 GeV. Panel (b): The scaled dimensionless jet 
transport parameters G/T® as a function of medium temperature T from the best fitting region of the panel (a). Panel (c): The single hadron 
suppression factors R.44(pr) with couples of (of /T?, do/Te) = (1.4, 3.6) (red solid curve) and (o2 /T2, ĝo /TẸ) = (00, 5.6) (blue dashed 


curve) compared with PHENIX [1, 2] data. 


medium temperature increase. The going-down T depen- 
dence of ¢/T® stimulates an attempt for a Gaussian assump- 
tion for the T dependence of ¢/T°. Here we assume the apex 
of the Gaussian distribution is seated at the critical tempera- 
ture, as shown in Eq. (11). This assumption for the Gaussian 
temperature dependence will be submitted into Eq. (13) with 
dp = O for the QGP phase only, and the Eq. (7) and (8) for 
jet energy loss. 


A. Fit R44 and v2 at RHIC 


In Eq. (11) for the assumption of the Gaussian temper- 
ature dependence, the introduced ĝo /T is still the scaled 
jet transport parameter in the initial time at the center of 
the medium, while o7,/T? is the squared Gaussian width. 
Starting with Au + Au collisions at Vsnn = 200 GeV, 
we choose o7./T? € [0.35,3.5] with bin size 0.35 and 


ĝo/TÈ € [0.4, 5.2] with bin size 0.4 and get 130 couples of 
(o2,/T?, Go/T@) for Eq. (11) and (13) with ân = 0. We firstly 
get the R4,4(pr) in 0 - 5% collisions and compare with data 
[1, 2] to get a 2-dimensional contour plot for the x?/d.o.f 
fitting, as shown in Fig. 9 (a). According to the obtained 
best fitting region for (o2 /T2, Go/T3), we give G/T? as a 
function of T in Fig. 9 (b). The red solid curve is for a going- 
down g/T® dependence with (o2 /T2, Go/Tj) = (1.4,3.6) 
while the blue dashed curve is for a constant dependence 
with (07,/T?, do/T3) = (00,5.6). Such two dependence 
forms give almost the same R44(pr), as shown in Fig. 9 
(c). Raa(pr) doesn’t “care” whether g/T? is of Gaussian 
temperature dependence or not at RHIC. 

Shown in Fig. 10 (a) and (b) are the y?/d.o.f analy- 
ses of single hadron R4,(pr) and elliptic flow v2(pr) as 
a function of o7,/T? and go/T} from fitting to experimen- 
tal data [1-3] in 20 - 30% Au + Au collisions at \/snn = 
200 GeV, respectively. The global x?/d.o. f fitting results for 
both R44(pr) and v2(pr) are shown in Fig. 10 (c). Al- 
though y?/d.o.f(R.aa) performs inactively for the tempera- 
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Fig. 11. (Gaussian-dependence) In 20 - 30% Au + Au collisions at ,/snn = 200 Gey, the scaled dimensionless jet transport parameters 
g/T® as a function of medium temperature T from the best fitting region of global x? fits of Fig. 10 (c) are shown in panel (a). The 
single hadron suppression factors R.44(pr) and elliptic flow parameter v2(pr) are shown in panel (b) and (c), respectively, with couples of 
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ture dependence in 20 - 30% as well as in 0 - 5% centrality, 
x? /d.o.f (v2) expresses a great favor in the small Gaussian 
width which gives G/T? going down more rapidly with T. 
Consequently, the global x?/d.o. f fits for both R.4a(pr) and 
v2(pr) give the explicit constraint for the introduced param- 
eters (o7,/T?, ĝo /T8). 

With the best global fitting values for (o7./T?, Go/Tj), we 
show in Fig. 11 (a) for the Gaussian T dependence of ¢/T®. 
Choosing the constant dependence with (02 /T2, Go/T}) = 
(00, 5.6) (blue dashed curve) and a Gaussian T dependence 
with (o%,/T?, Go/T@) = (1.4, 3.6) (red solid curve) for G/T?, 
we again get the almost same R44(pr) as shown in Fig. 
11 (b), and v2ə(pr) with difference less than 5% in Fig. 11 
(c). Compared to the case of the linearly-decreasing T de- 
pendence, the two v2(pr) in Fig. 11 (c) are closer to each 
other because the difference between the Gaussian T depen- 
dence and the constant dependence is smaller, which leads to 
an almost invisible change in v2(pr). 

The use of a Gaussian form for G/T was intended to pro- 
vide a more flexible temperature dependence by narrowing 


(1.4, 3.6) (red solid curve) and (co, 5.6) (blue dashed curve) compared with PHENIX [1-3] data. 


the Gaussian width compared to the linear form. Neverthe- 
less, although v2 data favor a higher Gaussian peak, while 
R44 fitting imposes constraints on it. The global x? fitting 
results for the Gaussian temperature dependence hypothesis 
presented in Fig. 11(a) did not manifest a steeper decline 
with increasing temperature. When comparing the red curve 
in Fig. 11(a) with that in Fig. 3(a), it is apparent that the lin- 
ear temperature dependence results in a higher q./T? at the 
critical temperature. And v2 is more sensitive to energy loss 
near the critical temperature as stated previously. Therefore, 
the performance of the Gaussian shape is not much better than 
that of the linear shape. 


B. Fit Raa and v2 at the LHC 


Make the same process going through the case for Pb+ Pb 
collisions at ,/snn = 2.76 TeV. For the most central col- 
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Fig. 12. (Gaussian-dependence) Panel (a): The x? /d.o. f analyses for single hadron R44 (pr) as a function of G/T? and ĝo/T9 from fitting 
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lisions, we choose o7./T? € [0.35,3.5] with bin size 0.35 
and ĝo /T9 € [0.2, 2.4] with bin size 0.2 and get 120 couples 
of (07,/T?, Go/T§) to get Raa(pr) and the 2-dimensional 
x?/d.o.f fitting contour plot, as shown in Fig. 12 (a). With 
the constant dependence with (02 /T2, ĝo/T8) = (co, 3.6) 
and the Gaussian T dependence with (07./T?,4o/T3) = 
(1.75, 1.6), we get the same R4,(pr) to fit data well, as 


shown in Fig. 12 (b). 


For 20 - 30% collisions, we choose o7/T? € [0.35, 7.0] 
and ĝo/T € [0.2,4.4], and separately get 440 groups of 
Raa(pr) and v2(pr) to perform the y?/d.o.f fitting, as 
shown in Fig. 13 (a) and (b), respectively. The global fits 
for both R4 (pr) and v2(pr) are shown in Fig. 13 (c) where 
the constraints for the introduced parameters (07,/T?, Go/T@) 
are obtained. With the constraints, the Gaussian T' depen- 
dence of å/T? is shown in Fig. 14 (a). Choosing a Gaussian 
T dependence with (o2 /T2, Go/T3) = (1.05, 1.2) (red solid 
curve) and the constant dependence with (07./T?, Go/T3) = 
(co, 4.1) (blue dashed curve) for ¢/T?, we again get almost 


the same R44(pr) shown in Fig. 14 (b), and different v2(pr) 
in Fig. 14 (c). Compared to the constant case, the going-down 
Gaussian T dependence of ĝ/T? gives v2(pr) an enhance- 
ment by about 10%. 


Both at RHIC and the LHC, numerical results for simulta- 
neously fitting to R44(pr) and v2(pr) show that the Gaus- 
sian T dependence of ¢/T® is smoothly going down with T 
and similar to the linearly-decreasing T dependence of ¢/T®. 
Compared with the constant ¢/T®, the going-down T depen- 
dence of @/T® enhances the hadron azimuthal anisotropy by 
about 5-10% to improve v2(pr) to fit data. 


Thus far, we have demonstrated the constraining power 
of the experimental data on the three temperature-dependent 
forms of ¢/T?. To more clearly distinguish the separate con- 
straining effects of R44 and v2, we have listed the best-fit 
parameters and the corresponding minimum y?/d.o.f values 
for each scenario in Table 1. These values correspond to the 
results shown in Figs. 2, 5,10, and 13. Note that for the con- 
stant form of ¢/ T?, the və data were not utilized. 
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Fig. 14. (Gaussian-dependence) In 20 - 30% Pb + Pb collisions at ,/snn = 2.76 TeV, the scaled dimensionless jet transport parameters 
g/T® as a function of medium temperature T from the best fitting region of global x? fits of Fig. 13 (c) are shown in the panel (a). 
The single hadron suppression factors Raa(pr) and elliptic flow v2(pr) are shown in panels (b) and (c), respectively, with couples of 
(a2 /T?, ĝo /TẸ) = (1.05, 1.2) (red solid curve) and (co, 4.1) (blue dashed curve) compared with experimental data [4, 5, 9, 10]. 


Constant Linear T dependence Gaussian T dependence 
Ge]Te = 90/15 X [4.0] Ge/Te, GoTo) x7/4.0.f| (07 /T?, Go/T3) x7 /d.0F 

20-30% Au + Au collisions at 200 GeV 

Raa (5.8,5.8) 0.17 (6.9,1.1) 0.17 (2.8,4.4) 0.48 

v2 = - (8.7,0.5) 2.39 (0.4,2.4) 2.98 

(Raa + v2) - a (7.5,0.2) 1.51 (0.7,2.4) 2.45 
20-30% Pb + Pb collisions at 2.76 TeV 

Raa (4.0,4.0) 0.59 (5.0,3.6) 0.54 (4.8,3.2) 0.56 

v2 = - (5.4,0.2) 1.52 (2.1,2.4) 1.23 

(Raa + v2) — = (5.1,0.2) 0.94 (1.4,1.6) 1.23 


TABLE 1. Optimal parameter and corresponding  /d.o.f for different data sets in different G/T? form. Note that for the constant form of 


G/T", the v2 data were not utilized. 


V. LINEAR TEMPERATURE DEPENDENCE OF ¢/T? IN 
QGP AND HADRON PHASES 


In this section, for a total study of the linear temperature 
dependence of air, we use Eqs. (10) (13) and (14) to in- 
clude contributions from both QGP and hadron phases in 20 
- 30% A + A collisions. 


A. Best fits for R44 and v2 due to energy loss in QGP and 
hadron phases 


Let us begin with Au + Au collisions in 20 - 30% at 
/8nN = 200 GeV. Shown in Fig. 15 are hadron suppression 
factors Ra a(pr) (panel (a)) and elliptic flow ve(pr) (panel 
(b)) with jet energy loss of both QGP and hadronic phases 
(black solid curves) and of only QGP phase (red solid curves 
from Fig. 3) for 20 - 30% Au-+ Au collisions at ,/snn = 200 
GeV, respectively. With the same go /T = 0.2 shown in Fig. 
3, we have to decrease the value of g./T? from 6.9 to 5.7 
due to the included hadronic phase to get the same R.44(pr). 
The reduction of 20% for G/T® is consistent with a study in 


Ref. [23]. The jet energy loss in the hadron phase gives an 
enhancement by 10% for elliptic flow v2(pr), which means 
the jet energy loss of the hadronic phase has an important and 
non-negligible contribution to v2(pr). 

As shown in Fig. 16 is for the case of 20 - 30% Pb + 
Pb collisions at ,/snn = 2.76 TeV. With the consistence 
Raa(pr) in Fig. 16 (a), the ve(pr) (black solid curve) 
with (@./T?,Go/T3) = (4.5,0.0) of QGP phase and qj, of 
hadronic phase has an additional enhancement by 5% com- 
pared with (@-/T?, Go/T$) = (5.4,0.0) of only QGP phase 
(red solid curve from Fig. 6) in Fig. 16 (b). The enhancement 
is less than that at RHIC because the fraction of hadron phase 
contribution to jet energy loss at the LHC is less than that at 
RHIC, as shown in Fig. 17 (a) and (d). 

It is worth noting that the description of pr dependence 
of v2 still deserves further improvement. This work only 
considers the T dependence of @/T?. Taking into account 
the dependence of g/T? on parton energy Æ, which is more 
directly related to the pr distribution of the jet energy loss 
[48, 50, 88], is likely to aid in improving the pp dependence 
of v2. Moreover, considering the dependence of å/T? on the 
colliding energy \/snn could also contribute to a more ac- 
curate description [89]. Additionally, elastic energy loss and 
jet-induced medium responses have a significant impact on 
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Fig. 15. Single hadron suppression factors R4 4(pr) (panel (a)) and elliptic flow v2(pr) (panel (b)) with jet energy loss of QGP + Hadron 


phases (black solid curves with (q./T?,Go/T@) = 


(5.7, 0.0) and gp) and of only QGP phase (red solid curves with (G./T?,Go/T3) = 


(6.9, 0.0) from Fig. 3 ) for 20 - 30% Au + Au collisions at \/snn = 200 GeV. 
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Fig. 16. Single hadron suppression factors R44 (pr) (panel (a)) and elliptic flow v2 (pr) (panel (b)) with jet energy loss of both QGP + Hadron 


phases (black solid curves with (ĝe /T3, ĝo /TÈ) = 


(4.5,0.0) and ĝa) and of only QGP phase (red solid curves with (ĝe /T3, ĝo/T9) = 


(5.4, 0.0) from Fig. 6) for 20 - 30% Pb + Pb collisions at \/snn = 2.76 TeV. 


hadron production at the intermediate transverse momentum, 
which could enhance the və at this pr region [90, 91]. We 
would incorporate these considerations into our model in the 
future to further improve its descriptive power for experimen- 
tal data. 


B. Jet energy loss distribution in QGP and hadron phases 


Similarly to Fig. 7 and 8, we show in Fig. 17 the com- 
parisons of the average jet energy loss distribution between 
QGP + Hadron phases (black solid curves) and only QGP 
phase (red solid curves) in noncentral A + A collisions at 
RHIC (left panels) and the LHC (right panels), respectively. 
The black dot-dashed curves are for hadron phase contribu- 
tion in the case of QGP + Hadron phases. From top to bottom 
are the average accumulative energy loss, differential energy 
loss, and the energy loss asymmetry, respectively. 

The hadron phase contribution to the total energy loss of 
QGP + Hadron phases is about 17% at RHIC in Fig. 17 (a) 
while about 14% at the LHC in Fig. 17 (d). Because of the 
first order phase transition in the current model study shown 


in Fig. 7 (c) and (f), the hadron phase contribution happens 
mainly in the T, nearby shown in Fig. 17 (b) and (e), which 
strengthens the azimuthal anisotropy of the system and then 
enhances the elliptic flow parameter. This is similar to the 
peak of the energy loss rate pushed to the T, due to the going- 
down linear T dependence of åĝ/ T? in Fig. 7 (b) and (e). Such 
enhancements of the azimuthal anisotropy are also exhibited 
by the energy loss asymmetry shown in Fig. 17 (c) and (f). 


C. Temperature dependence of ĝâ/T? in QGP and hadron 
phases 


Shown in Fig. 18 is the ¢/T? of the hadronic phase (dot- 
dashed line) and QGP phase (solid lines) as a function of 
medium temperature T. The hadronic phase ¢/T® is given 
by Eq. (14). The values of ¢/T? of QGP phase with the 
linear temperature dependence are given by Eq. (10) with 
(G-/T?, Go/Tj}) = (5.7, 0.0) at RHIC (denoted in red curve) 
and (4.5,0.0) at the LHC (denoted in blue curve), respec- 
tively. The contributions to ¢/T? of QGP and hadron phases 
combined via Eq. (13) are applied for R44(pr) and v2(pr) 
at RHIC/LHC in Fig. 15 and Fig. 16. As comparisons, the 


5.00 
200 GeV, Au+Au -->7°, 20-30% 
4.00 H7 
= E**=10 GeV 
% 3.00 
g 
a 2.00 
3 
1.00 
0.00 
= only QGP phase (b) 
E (G./T2, Go/T3)=(6.9, 0.0) 
> 2.00; QGP + Hadron phases 
2 (Q/T2, Qo/T2)=(5.7, 0.0) and ĝ, 
G P 
K 1.00 hadron phase contribution 
< 
g 


a 
[e] 
(=) 


0.10 


A (AE) 


0.05 


q(fm/c) 


14 


= 
he 
fo) 


2.76 TeV, Pb+Pb -->h*, 20-30%] (d) 


EJet=100 GeV 


1o 
ro) 


(AE) (GeV) 
a 
[=] 


only QGP phase (e) 
(G/T, Âo/T3)=(5-4, 0.0) 

QGP + Hadron phases 
(G/T2, Qo/T2)=(4.5, 0.0) and ĝ, 
hadron phase contribution 


(dAE/dt) (GeV/fm) 


t(fm/c) 


Fig. 17. Comparisons of the average jet energy loss distribution between QGP + Hadron phases (black solid curves) and only QGP phase 
(red solid curves from Fig. 7 and 8) in noncentral A + A collisions at RHIC (left panels) and the LHC (right panels), respectively. The black 
dot-dashed curves are for hadron phase contribution in the case of QGP + Hadron phases. From top to bottom are the average accumulative 
energy loss, differential energy loss, and energy loss asymmetry, respectively. 


representative samples of å/T? extracted from single hadron, 
dihadron, and y-hadron production at RHIC and the LHC en- 
ergies with information field (IF) Bayesian analysis [51] are 
also shown in Fig. 18, which are denoted by grey curves. Ref. 
[51] do not consider the jet energy loss in the hadronic phase 
with a pseduocritical temperature T, = 0.165 GeV. When 
considering the experimental data covering a wider tempera- 
ture range, the g/T® still points out a large value at the crit- 
ical temperature. Our numerical results for the temperature- 
dependent ĝ/T? are consistent with that obtained by the IF- 
Bayesian method. After adding the hadron phase contribu- 
tion to jet energy loss, one should decrease the QGP phase 
contribution so as to get the total energy loss equal to that for 
the case of only the QGP phase. The decreased QGP-phase 
energy loss makes v2(pr) smaller while the added hadron- 
phase energy loss makes v2(pr) larger. Numerical results 
show that the competition between them gives v2(pr) larger 
than the case of only the QGP phase due to the energy loss 
contribution of the hadronic phase concentrating near the crit- 
ical temperature. No matter the added hadron phase or the 
going-down linear T dependence of ¢/T® in the QGP phase, 
as shown in Fig. 18, the jet is made to lose much more energy 
near the critical temperature, which results in a larger v2(pr) 
for the large pr to fit data better. Such a numerical result 
of the stronger jet quenching in the near-Te region is consis- 
tent with a previous theoretical study [33]. One may notice 
that we extract the g/T® at RHIC and the LHC with differ- 


ent parameter ranges. This is mainly because, in this work, 
we observed that extracting ĝâ/T? at RHIC and the LHC sep- 
arately could give a better x? result compared to the simulta- 
neous fit at both colliding energies. We only considered the 
dependence of ĝ on temperature in this work. As mentioned 
previously, if the dependencies of ĝ on parton energy E and 
\/SNN; as well as the jet-induced medium responses and the 
elastic energy loss, were all taken into account, the value of 
q could be constrained more accurately from R44 and v2 
at both RHIC and LHC energies simultaneously. Neverthe- 
less, the ¢/T? obtained in this work is also consistent with 
the ĝâ/T? extracted by the IF-Bayesian approach [51] and the 
JETSCAPE [48, 50], which constrained the ĝ/ T? from RHIC 
and the LHC simultaneously, as shown in Fig. 18. We hope 
that in the future, by updating the model and enriching the 
information on g, we can better simultaneously describe the 
Raa and v2 for both RHIC and the LHC energies. 


VI. SUMMARY 


In this paper, within a next-to-leading-order perturbative 
QCD model, the medium temperature dependence of jet en- 
ergy loss is studied via nuclear modification factor R4 (pr) 
and elliptic flow parameter vo(pr) of large transverse mo- 
mentum hadrons. Due to jet quenching, medium-modified 
fragmentation functions based on the higher-twist energy loss 
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Fig. 18. The scaled jet transport coefficient @/T? as a function of 
medium temperature T. T < 0.17 GeV for the hadronic phase (dot- 
dashed curve), and T > 0.17 GeV for the QGP phase (solid curves). 
The g/T® of QGP phase for Au + Au collisions at J/snn = 200 
GeV is denoted in red and for Pb + Pb collisions at ,/snn = 2.76 
TeV is denoted in blue. As comparisons, the ĝ/ T? posterior samples 
from single hadron, dihadron, and y-hadron production at RHIC and 
the LHC energies with IF-Bayesian analysis [51] are also shown in 
grey curves. 


formalism are used in the numerical calculations. We assume 
the scaled jet transport coefficient å/T? depends on medium 
temperature by linear or Gaussian form, with which we have 
calculated the single hadron suppression factor R44 (pr) and 
the elliptic flow parameter v2(py) and compared them with 
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the LHC while go /T$ = 0.0-4.2 at both RHIC and the LHC, 
as shown in Fig. 3 (a) and Fig. 6 (a). The numerical results 
show that R44(pr) and v2(prp) are both more sensitive to the 
value of G/T? near critical temperature T, than near initial 
highest temperature To. Furthermore, the fitting results give 
a trend of the going-down q/T® dependent on the medium 
temperature, which is also supported by the Gaussian T de- 
pendence of ¢/T®? for only the QGP phase. 

Compared with the case of constant ¢/T®, the going-down 
T dependence of ¢/T? makes a hard parton jet to lose more 
energy near the critical temperature T, and therefore strength- 
ens the azimuthal anisotropy for large pr hadron produc- 
tions. As a result, the elliptic flow parameter v2(pr) for large 
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RHIC/LHC. Considering the first order phase transition from 
QGP to hadron and the hadron phase contribution to jet en- 
ergy loss, v2(pr) is again enhanced by 5-10% at RHIC/LHC. 
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